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Abstract. The paper presents a systematic derivation of macroscopic equations for 
freeway traffic flow from an Enskog-like kinetic approach. The resulting fluid-dynamic 
traffic equations for the spatial density, average velocity, and velocity variance of ve- 
hicles are compared to equations, which can be obtained from a microscopic force 
model of individual vehicle motion. Simulation results of the models are confronted 
with empirical traffic data. 

1 Introduction 

During the last five years, modelling and simulating traffic dynamics has found 
a large and rapidly growing interest in physics. This is due to 

1. similarities of traffic dynamics with fiows of gases, fluids, and granular media, 

2. instability phenomena and critical behavior of traffic (cf. Fig. 1), 

3. interesting applications of cellular automata and molecular dynamics simu- 
lation methods, 

4. the need of efficient traffic optimization methods in order to keep or increase 
the level of mobility. 

Usually, one distinguishes three levels of modelling: The microscopic level of de- 
scription delineates the dynamics of the single driver- vehicle units [1] , [2] , [3] , [4] , 
[5]. This allows to consider different vehicle characteristics and driving styles, so 
that many of these models aim at a high-fidelity description of traffic flow, e.g. 
[3], [5]. They are mostly used for detail studies (e.g. of on-ramp traffic, bottle- 
necks, effects of traffic optimization measures), but they consume an enourmous 
amount of CPU time due the the large number of variables involved. An alter- 
native approach are cellular automata, which allow to simulate a minimal model 
of traffic dynamics faster than real-time [6] , [7] , [8] , [9] . 

Computational efficiency can also be reached by macroscopic traffic models, 
but at a higher degree of accuracy [10], [11], [12], [13], [5]. Macroscopic traffic 
models consist of equations for a few aggregate quantities like the spatial den- 
sity p, the average velocity V , and (in some cases) additional velocity moments. 
These equations are similar to fluid-dynamic equations, but some fundamental 
differences with respect to the dynamics of ordinary fluids have recently been 
recognized [5], [13]. For congested conditions, their detailled form is not at all 
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Fig. 1. Temporal evolution of the average velocity V{r,t) at subsequent cross-sections 
of the Dutch highway A9 from Haarlem to Amsterdam at October 14, 1994 (five minute 
averages of single vehicle data) [14], [5]. The prescribed speed limit is 120km/h. We 
observe a breakdown of velocity during the rush hours between 7:30 am and 9:30 am 
due to the overloading of the highway at r = ro := 41.8 km (■ ■ ■). At the subsequent 

cross-sections the traffic situation recovers ( : r = ro + 1km; : r = ro + 2.2 km; 

— : r = ro + 4.2 km). Nevertheless, the amplitudes of the small velocity fluctuations at 
ro become larger and larger, leading to so-called stop-and-go waves. 

obvious. Therefore, it has been suggested to derive the macroscopic traffic equa- 
tions from a kinetic, i.e. mesoscopic level of description, which delineates the 
spatio-temporal evolution of the velocity distribution [5], [13], [15], [16], [17], [18]. 

2 Enskog-like Kinetic Traffic Model 

In the following, we define the coarse-grained phase-space density p{r,v,t) of 
vehicles per lane with velocity v at place r and time t by 

p{r,v,t) := (2z\r)^(2z^^) E J dr' J dv' S{r' - r^{t)) S{v' - v^{t)) . (1) 

ra{t) is the location and Va{t) the velocity of vehicle a at time t. We do not 
distinguish different lanes, here, although this is possible [5], [19]. Instead, we 
treat the overall cross section of an n-lane freeway in an effective way [5] , [20] . 



From Microscopic to Macroscopic Traffic Models 3 



Let us assume an acceleration equation of the form 

'^=h{Va)+ ^ fa(3{roc,Vcc-,ri3,V(3)+ia{t)-, (2) 



dt 



where the function 



r , ^ Vq- Va{t) , , 

fo{Va) ■■= — (3) 



describes an adaptation of the actual velocity Va{t) to the desired velocity Vq 
within a (possibly density-dependent) relaxation time r. The second term in 
(2) delineates the effect of interactions with vehicles (3, and ^a{t) reflects ve- 
locity fluctuations due to imperfect driving. Wc will assume {^a{t)^(j{t')) = 
2D 5aj3 5{t—t'), where the diffusion function D is density- and velocity-dependent 
[5], [13], [26]. For reasons of simplicity, the desired velocity Vq and the relaxation 
time T were taken identical for all vehicles, but it is also possible to generalize 
this model [16], [17], [18]. 

Prom (1) and (2), the following dynamical equation for the phase-space den- 
sity can be derived: 



The fluctuation term gives a contribution 

'dp\ d\pD) 



dt J jj dv"^ 



In addition, we have used the abbreviation 



(5) 



with the average interaction force 



'''' lm«} (6) 



r+Ar v+Av 



/int(r, V, t) := ^~2rAv ^ J J ^""^ ~ """^^^^ ~ ' 

"^^r—Ar v—Av 

The interaction term (6) reflects deceleration processes. In analogy to the En- 
skog theory of dense gases [21] and granular media [22], [23], [24], but with an 
interaction law typical for vehicles [5], [13], it is approximated by 

={l-p)x{r + l,mv) (8) 
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with the Boltzmann-like interaction function 

B{v) = J dw {w — v) p{r,w,t)p{r + s,v,t) 

w>v 

— J dw {v — w)p{r,v,t)p{r + s,w,t) . (9) 

v>w 

According to this, the phase-space density p{r, v, t) increases due to deceleration 
of vehicles with velocities w > v, which cannot overtake vehicles with velocity 
V. The density-dependent probability of immediate overtaking is represented by 
p. A decrease of the phase-space density p{r, v, t) is caused by interactions of 
vehicles with velocity v with slower vehicles driving with velocities w < v. The 
corresponding interaction rates are proportional to the relative velocity |?; — w| 
and to the phase space densities of both interacting vehicles. By s{V) = lo + l{V) 
(w vehicle length -I- safe distance) it is taken into account that the distance 
of interacting vehicles is given by their velocity-dependent space requirements. 
These cause an increase of the interaction rate, which is described by the pair 
correlation function xi"!") = [1 ~ pI'^j^)*]"^ at the 'interaction point' r + I. A 
more detailled discussion of the above kinetic traffic model is presented elsewhere 
[5], [13]. 

Now, we will focus on the the macroscopic equations for the spatial density 
p{r,t)= jdvp{r,v,t), (10) 

the average velocity 

''<^'" = /<'-^' 

and the velocity variance 

e{r,t)^ Jdv[v-V{r,t)rP^^. (12) 

These are obtained by multiplying the kinetic equation with f and integrating 
with respect to v. In order to obtain a closed system of equations, we assume 
that the velocity distribution P{v; r, t) has a Gaussian form, i.e. 

d(r V t) e-[^-^('''*)l'/P«('''*)l 
P(v. r, t) := PSLl^ = t . (13) 

According to empirical data, this approximation is well justified (cf. Figs. 2 and 
3). 

After some straightforward calculations, the following equations are obtained: 
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Fig. 2. Comparison of empirical velocity distributions at diflterent densities ( — ) 
with frequency polygons of grouped Gaussian velocity distributions with the same 
mean value and variance ( — ) [25]. A significant deviation of the empirical relations 
from the respective discrete Gaussian approximations is only found at a density of 
p = 40 vehiclcs/km lane, where the temporal averaging period of T = 2min may have 
been too long due to rapid stop-and-go waves. 



^^^^K^[p{V^+e)] = ^{Vo-V) + {l-p)x{r + l,t) JdvvB{v), (15) 



dt 



+ (1 - p)x(r + l,t) j dv v^B{v) . (16) 



Equations (14) to (16) are similar to the Euler equations of ordinary fluids. In 
particular, the density equation (14) agrees with the continuity equation, reflect- 
ing that the number of vehicles is conserved (on a circular road). However, equa- 
tions (15) and (16) contain some additional terms compared to the hydrodyamic 
equations for momentum density and energy density, which are essential for the 
instability of traffic flow. The respective first terms on the right-hand sides of 
(15) and (16) originate from the acceleration towards the desired velocity Vq and 
from velocity fluctuations. The respective last terms reflect interaction (deceler- 
ation) effects. In contrast to ordinary fluids, they do not vanish, since vehicular 
interactions do not conserve momentum and energy. 
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Fig. 3. Density-dependence of the skewness 7 of the velocity distribution (• : 1-minute 

data; O: respective mean values) [25]. The largo variation of the 1-minuto data at low 
densities is due to the small number of vehicles which pass a cross section during the 
time interval T = 1 min, whereas the large variation of their mean values at high den- 
sities comes from the few 1-minute data, over which could be averaged. The 1-minute 
data of the skewness scatter around the zero lino ( ) and mostly lie between —1 and 
1, so that it is negligible most of the time. 



To obtain the explicit form of the interaction terms, one has to carry out a 
number of lengthy calculations. Using the abbreviations 



p+{r,t):=p{r + s,t), V+{r,t) := V{r + s,t) , e+{r,t) := e{r + s,t) , (17) 
and introducing the Gaussian error function 



^{x) = / dy 



2-K 



(18) 



one finally finds 



j dvvBiv) = -PP+ I [{0 + 6+) + {V- V+f] $ {j== 



+ {v-v+){e + e+)- 



-{V-V+f/[2(0+B+)] 



(19) 
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and 



/ 



dvv^B{v) = -2pp+{0-0+) 



{o + o+y 



-{v-v+f/[2{e+e+)] 



+{V-V+) <P 



^2TT{e + 6+) 

+ {V + V+) I dvvB{v). (20) 



The macroscopic traffic equations (14) to (16) were; written as equations for 
fluxes with sink/source terms (the terms on the right-hand side). The flux rep- 
resentation is very advantageous, since many numerical integration algorithms 
have been developed for this type of partial differential equations. However, due 
to (19) and (20), the flux representation is non-local. This is caused by the flnite 
space requirements of cars, i.e. a driver reacts to another car already at a cer- 
tain distance. As a consequence, the non-local interaction terms imply viscosity 
effects, among other things. To see this, we expand them up to second order. 
Neglecting products of spatial derivatives, we get the continuity equation 



dp ^dp dV 
dt dr ^ dr 



(21) 



the velocity equation 



dV_ ydV_^ dp 
dt dr dr 



02 



dV 
dr 



de 

dr 



, 5V , d^V , d^0 Ve 



V 



dr^ 



Qj.2 



Q^2 



(22) 



and the variance equation 



de ^de _ ^ dp 
dt dr dr 



C2 



di 



9V 



dV_ 
dr 
d^V 



d0 
dr 



d2 



Qj.2 



d^0 2{0e - 0) 
^ dr'^ T 



(23) 



(which corresponds to the equation of heat conduction in ordinary fluids). Here, 
we have used the abbreviations 



ai = 


-[^ 


+ (l-p)xs(l 




a2 


= (1- 


P)xp(2s^- 


as = 


-[1 


+ (i-p)pxi] 




bi 


= -(1 


-p)x«(f + PX 


b2 = 


(1- 


p)xp{f\[^- 


P^Xw) ' 


b3 


= -(1 


-p)xpt ' 


C2 = 


-[2 


+ {l-p)Xps\0 




cs 


= 2(1- 


-P)XPS\/^, 


d2 = 


-(1 


-P)XPS^I > 




d3 


= (1- 


P)XPS^\/^, 



(24) 



and 



V, = Vo-t{1-p)px0, 0, = tD. 



(25) 



8 Dirk Helbing 



Note that ci = and di = 0, which is a consequence of the assumed interaction 
law of vehicles. It is one of the advantages of a kinetic derivation of macroscopic 
traffic equations, that the above functions can be analytically calculated. For 
example, we have obtained an expression for the equilibrium velocity Vc- Ac- 
cording to (25), it is given by the desired velocity Vq, diminished by a term due 
to decelerating interactions. The latter is proportional to the vehicle density and 
to the velocity variance, which is very plausible. The function dV /dp := —pai 
can be interpreted as the partial derivative of the "traffic pressure" V with re- 
spect to density. The quantity i] := pb2 has the meaning of a viscosity, which 
smoothes out sudden spatial changes of the velocity profile V{r, t). Both quanti- 
ties are non- negative and diverge at maximum density Pmax := l/^o, as it should 
be for reasons of consistency [5], [13]. Previous macroscopic traffic models did 
not describe these important facts correctly, since they have neglected the terms 
in (24) which explicitly contain I or s. Therefore, they are not valid for high 
vehicle densities, i.e. for congested conditions. Finally, note that it is possible to 
calculate Navier-Stokes corrections of the coefficients aj, hi, Ci, and di [26]. 




Fig. 4. Instability diagram for the Euler-like macroscopic traffic equations, including 
the dynamic variance equation [26] . 

According to our approximations, equations (21) to (23) are valid for small 
gradients of p, V, and 9. Therefore, they allow to investigate the evolution of 
small disturbances of the (stationary and spatially homogeneous) equilibrium 
solution. Figure 4 depicts the result of a linear instability analysis, showing that 
traffic flow is only stable at small and extreme densities as well as large wave 
numbers jfcj (i.e. small wave lengths £ = 27r/|fc|). This is in agreement with 
empirical findings. 

The instability diagram is obtained by 
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1. assuming a small periodic perturbation Sg{r, t) — go exp[ifcr + (A + iLu)t\ of 
the macroscopic traffic quantities g G {p, V, 9} relative to the stationary and 
spatially homogeneous equilibrium solution gdp) (go being the amplitude, k 
the wave number, A the growth rate, and uj the frequency of the perturba- 
tion), 

2. inserting g{r, t) — g^ + 6g{r, t) into the macroscopic traffic equations, 

3. neglecting quadratic terms in the small perturbations 5g{r,t) <^ g^, 

4. determining the three complex eigenvalues A = A + icj of the linearized 
equations in dependence of p and k. 

An explicit example for this procedure is discussed in [5], [13]. Equilibrium traffic 
flow is unstable if at least one of the growth rates is positive, i.e. max A > 0. 
Therefore, the instability diagram shows max A(/c, p) if this is greater than zero, 
otherwise 0. Figure 5 depicts a simulation result which demonstrates emerging 
stop-and-go traffic. 




5 



Fig. 5. Above a certain density, traffic flow is unstable. This gives rise to the develop- 
ment of stop-and-go waves. In the represented simulation, we applied periodic boundary 
conditions (which corresponds to a circular road). 

2.1 Non-linear Phenomena 

As a consequence of the inherent non-linearity of the macroscopic traffic equa- 
tions, they display the self-organization of a number of collective patterns of mo- 
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tion. This includes the formation of density clusters ('traffic jams'), anti-clusters, 
dipole layers, cascades of density clusters ('stop-and-go traffic'), and merging of 
density clusters. Moreover, one finds subcritical instabilities and non-linear wave 
selection phenomena [11], [12]. 



3 An Alternative Approach 

The Boltzmann-like formula (9) for vehicle interactions implicitly assumes, that 
deceleration maneuvers happen instantaneously. This approximation is only valid, 
if the average duration of deceleration maneuvers is considerably smaller than 
the time scale of the macroscopic traffic dynamics. However, one can also derive 
fluid-dynamic traffic equations without this assumption. We will illustrate this 
for a one-lane microscopic traffic model without possibilities of overtaking. 

3.1 A Concrete Microscopic Model 

Let us start with the social force model of vehicle dynamics, given by dva/dt = 
Va{t) and 

-TT = +/aln + l)('-o-'V,-/n + l-'V, + l)+^a(t)- (26) 

at T ^ ^ 

^ V ' " 

Acceleration Deceleration 

It is known that models of this kind are able to describe the emergence of stop- 
and-go traffic [5], [4] (cf. Figure 6). 

The advantage of the social force concept is, that it allows a very intuitive 
modelling of decision processes which are related to continuous changes in some 
(possibly abstract) space of behavioral alternatives [27]. According to this, the 
different motivations which influence an individual at the same time, are de- 
scribed by additive, force-like quantities. These generalized forces are, of course, 
no Newtonian forces. For example, they do not obey the law actio = reactio. The 
social force concept is well compatible with theoretical concepts from the social 
sciences and has been elaborated in detail [27], [5]. It has already been successful 
in describing various self-organization phenomena in pedestrian crowds [5], [28], 
but it was also applied to opinion formation processes [27], [29]. 

In the case of driver behavior, we have two conflicting motivations: On the one 
hand, the driver likes to accelerate towards his desired velocity Vq. On the other 
hand, he wants to keep a safe distance to the car in front. The latter is described 
by a repulsive deceleration force /«(»+!)• Effects /a/3 of interactions with other 
vehicles (3 ^ (a + 1) have been assumed to be negligible, here. However, they 
could easily be included in accordance with Eq. (2). 

As Fig. 7 shows, a good agreement with empirical data of driver-vehicle be- 
havior can be achieved with the following form of the repulsive interaction force: 

f ._ Vl{ra+i - rc) - , 

Ja{a+l) ■— r/Q(a-|-l) V^'J 
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Fig. 6. Representation of the vehicle trajectories of each 10th car on a circular road [5]. 
The slopes of the trajectories indicate the respective vehicle velocities, whereas their 
density reflects the spatial vehicle density. The simulation starts with a homogeneous 
traffic situation (i.e. all vehicles have initially the same distance to the car in front). 
In the course of time, density clusters emerge. These are often called 'phantom traffic 
jams', since they do not originate from any localized bottleneck. 



with 



f 

J a 



a(a+l) 



exp 



s(u„) 



R 



- Ua+l) . (28) 



Here, Qi^Av) is the Heaviside step function. 

If we would restrict the model to the first term of (27) (i.e. in the case r' ^ 
oo), we would arrive at the microsimulation model by Bando et al. [4]. V^{Ar) 
is the equilibrium velocity, which is a function of the distance Ar := r^+i — 
to the next car ahead. The additional term (28) takes into account that 



1 



Va+i is large 



drivers brake stronger, when the relative velocity Av := i 
or when the distance Ar to the car in front is small, 

2. the deceleration time r' is shorter than the acceleration time r, 

3. drivers begin to brake at a larger distance, if they drive fast. This is described 
by the velocity-dependence of the safe distance s{v). R is the range of the 
repulsive effect of a car. 

In most microsimulations, the relations Vl{Ar) and parameters t, t', R are 
specified individually (i.e. in an a-dependent way). 

It can be shown that the above force model is consistent in the limiting cases. 
For large distances or ~ Va+i, vehicle a approaches the distance-dependent 
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Fig. 7. Time-dependent velocity of a cax which follows another car in city traffic. The 
simulation results for the social force model treat the velocity of the vehicle ahead 
and its initial distance as given. We find a good agreement with the empirical fol- 
low-the-leader data. 



equilibrium velocity V^: 

dt ^ T 



(29) 



For small distances and Va > Va-\-i, it decelerates to the velocity Wa+i of the car 
in front: 

dVa Va+l{t) - Va{t) 



dt r'e['""+i~''»~*(''«)l/-^ 
With decreasing distance it brakes stronger. 



(30) 



3.2 Derivation of Macroscopic Traffic Equations 

In the following, we write the acceleration relation in the form 
dv V'(Ar) - V 

- = f{Ar, V, w) + m := ' + fi^r, v, w) + (31) 

at T 

with the abbreviations v := Va, w := Va+i, Ar := Va+i — Va, S, ■= ^a, and 



f '■— /a(Q+l) — /a(a+l) + 



(32) 
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Relation (31) is now inserted into the kinetic equation 

which is again a direct consequence of definition (1). Note that the interaction 
effects were absorbed into the function /, here. Next, we multiply this equation 
with v'^ and P'(w, Z\r|w, r, i), which denotes the probability that, given a car 
with velocity v is located at place r, the car in front drives with velocity w at a 
distance Ar. Finally, the resulting equation is integrated with respect to w and 
Ar. This gives the macroscopic equations 

dp d{pV)_ _ , . 

^ + |;[p(y2 + 9)] = P-{V: -V)+ pT^ , (35) 
+ Q)\ + ^^\p{y^ + 3V^0)] = ^(K*^ + rD - - 0) + pT^ (36) 

with 

jrfe(r,i) ■.= k j dAr J dv j dw v^-^f'{Ar, v, w)P'{Ar, w\r, v, t) ^^^^^'^^ 
= dAr J j dwv''~^ex.p^- ^^ i?^^^^ ) 



and 



xP'{Ar,w\r,v,t)P{v;r,t) (37) 



V* := j dAr j dv j dwV^{Ar)P'{Ar,w\r,v,t)P{v;r,t) . (38) 



To get (37), we made use of partial integration. 

In the following, we will apply the factorization approximation 

P'(Ar,w\r,v,t) P^{Ar;r,t)P{w;r + Ar,t) , (39) 

which is even exact, if the distributions of the velocities w and the headways 
Ar are statistically independent, and independent of v. Furthermore, we assume 
that the headway distribution P^,{Ar;r,t) is a function of the density p and 
average velocity V at a certain place r + Sr between r and r + Ar: 



P^Ar; r, t) = P^ (zir; p{r + 5r, t), V{r + Sr, t)) 



(40) 
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Then, we can expand (37) and (40) in the small quantities 5v := v — V and Sr, 
respectively. In this way, we obtain 

^/ X /", /", k ^v-w f Ar- siV) \ ( ds 5v 
W)^-kjdArjd. j c^«;.-^-^exp(^ R^){dVR 



+ 7 



dsY 
dV^ ' \dv) 



i?2 



P* {Ar; r, t)P{w; r + Ar, t)P{v; r, t) 



(41) 



and 



P. {Ar; r, t) ^ P. {Ar; p{r, ^) - ^(M)) + ^ + ^) 



'dV \~d^ ''^lh^~ 



(42) 



where wc again neglected products of partial derivatives dg/dr. After carrying 
out the integrations over w, w, and Z\r, the resulting macroscopic traffic equations 
can again be written in the form of Eqs. (21) to (23). However, the coefficients a^, 
hi, Ci, and rf,; arc different, since we did not apply the approximation of sudden 
deceleration maneuvers. The problem of this method is, that it does not provide 
the functional form of the headway distribution P^{Ar; p,V), which is needed 
for the explicit evaluation of the coefficients. Nevertheless, the use of the above 
results will be presented by a simple example. 



3.3 Relation between Bando's microscopic and Payne's macroscopic 
traffic model 

The microsimulation model by Bando et al. [4] is obtained by neglcicting the 
fluctuation term and the second term in (27), i.e. by setting D := and /' := 0. 
In order to calculate the corresponding macroscopic traffic equations, wc make 
a very simple assumption, here, namely that the headways Ar are given by the 
inverse of the density: 



P4Ar;r,t) := S [ Ar 



1 



p{r + dr, t) 



(43) 



Inserting this into the above equations, we finally arrive at the continuity equa- 



the velocity equation 

1 d{pe) 

p dr 
1 d{pe) 
p dr 



dt dr 



vn- 



-%{p) 



V 



1 dV* 

Tp^ dAr 



dp 
dr 



d'^p Sr"^ 



IdVofdp^ d^pSr^\ 
"^-^r^lTr'^+d^-)' ^''^ 



From Microscopic to Macroscopic Traffic Models 15 



and the variance equation 



d6_ ^de 
dt dr 



-29 



where 



:= v: 



dV_ 
dr 

1 
P 



T 



(46) 



(47) 



Close to the equilibrium solution, the variance equation can be neglected due to 
^ « 0. The instability condition of the remaining equations (44) and (45) reads 



dp 



! (5r 
> — 

T 



(cf. [5], [13]). This is only compatible with the instability condition 



dV* ! \_ 
dAr 2t 



or 



dV,. 



dp 



> 



2t 



of the Bando model [4] , if we choose 



6r = 



2p 



Ar 
~2 



(48) 



(49) 



(50) 



which is very plausible. In this case, the macroscopic equations (44) and (45) 
agree with the traffic model by Payne [10], but they contain the additional 
term ySr"^ j (2T)\{dVc I d p)d'^ p I dr"^ ^ which describes a smoothing of sudden spatial 
changes in density and velocity. However, as soon as the approximation « 
becomes invalid, Payne's model does not anymore reflect the traflic dynamics 
according to Bando 's model. 



4 Summary and Conclusions 

We have presented microscopic and macroscopic traffic flow models for freeways, 
which were successfully confronted with empirical data (cf. also [5], [13], [14]). 
Moreover, it has been shown, how macroscopic traffic models can be systemat- 
ically derived from the equations of motion for single vehicles. This is essential 
for increasing the speed of traffic simulations. Apart from the kinetic approach 
to this problem, which based on the assumption of sudden deceleration maneu- 
vers, an alternative method has been proposed, which presupposes a suitable 
approximation of the headway distribution function. The resulting macroscopic 
traffic equations are related to the hydrodynamic equations of ordinary fluids, 
but they contain a number of additional terms for three reasons: 

1. Vehicles accelerate to a certain desired velocity. 

2. A finite equilibrium variance of vehicle velocities is caused by imperfect driv- 
ing. 
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3. Vehicle interactions are anisotropic and do not conserve energy or momen- 
tum. 

The additional terms are responsible for certain instabilities of traffic flow, caus- 
ing 'phantom traffic jams' or 'stop-and-go trafSc'. They are also the origin of 
viscosity effects and of the divergence of 'traffic pressure' at maximum vehicle 
density. Here, it is essential that vehicular space requirements are taken into 
account [5], [13], [26]. Otherwise, the macroscopic traffic model would neglect 
certain characteristic terms, which would limit its validity to non-congested traf- 
fic situations. 

For the purpose of computer simulations, it is advantageous to have the 
macroscopic traffic equations in flux representation. This has been analytically 
derived, but it contains the Gaussian error function. In contrast to previous 
results [26], the corresponding equations are not restricted to cases of small 
gradients. 
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